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Abstract 

The effective action for QED in curved spacetime includes equivalence principle violating 
interactions between the electromagnetic field and the spacetime curvature. These interac- 
tions admit the possibility of superluminal yet causal photon propagation in gravitational 
fields. In this paper, we extend our analysis of photon propagation in gravitational back- 
grounds to the Kerr spacetime describing a rotating black hole. The results support two 
general theorems - a polarisation sum rule and a 'horizon theorem'. The implications for 
the stationary limit surface bounding the ergosphere are also discussed. 
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1. Introduction 

The possibility, originally discovered by Drummond and Hathrell in 1980 [1], that pho- 
tons propagating in curved spacetime may travel with speeds exceeding the usual speed 
of light has been the subject of renewed interest in the last two years[2-5]*. The phe- 
nomenon arises because vacuum polarisation in QED induces interactions between the 
electromagnetic field and spacetime curvature. Such interactions violate the strong prin- 
ciple of equivalence and allow the possibility of spacelike photon propagation without 
necessarily implying any violation of causality. 

In their original paper [1], Drummond and Hathrell studied photon propagation in the 
Schwarzschild spacetime, together with other examples including gravitational wave and 
Robertson- Walker backgrounds. In our previous paper [2], we extended this analysis to 
the Reissner-Nordstrom geometry describing a charged black hole. Here, we complete our 
survey of photon propagation in black hole spacetimes by considering the Kerr geometry 
describing a rotating black hole. 

Our study is motivated by the hope that by examining special cases with a particularly 
rich structure it may be possible to uncover general properties of photon propagation in 
gravitational fields. That hope is borne out by the results presented here. In particular, 
it had been observed previously that the light cone for radially directed photons in both 
the Schwarzschild and Reissner-Nordstrom geometries remains unperturbed, i.e. at k 2 = 0. 
For the Kerr spacetime, we find that this is no longer true and photons travelling on radial 
trajectories may indeed have velocities differing from 1, either greater or smaller than the 
usual velocity of light depending on the polarisation. However, we also find that these 
corrections always vanish at the event horizon itself. The light cone always remains k 2 = 
at the horizon. 

This result, together with a further observation that in all cases the corrections to 
the photon velocity are equal and opposite for the two transverse polarisations, was the 
motivation for ref.[12] where they are formalised into two theorems - a polarisation sum 
rule and a 'horizon theorem' (see section 5). This latter paper also contains some related 
observations about electromagnetic birefringence and the role of the conformal anomaly in 
photon propagation. 

2. Photon Propagation in Curved Spacetime 

In this paper, we consider the nature of photon propagation implied by the action 

T = Jdxy/^(-^ v F^+^(aRF^F^+bR^F^F\+cR^ (2.1) 

where the mass scale m and constants a, b and c are regarded as free parameters. This is 
the simplest action which incorporates explicit equivalence principle violating interactions 

* For earlier papers related to ref. [1] see refs.[6,7]. Similar phenomena in non-gravitational backgrounds 
are discussed in refs. [8-11,3]. 
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and leads to possible modifications of the light cone. (Other more complicated actions are 
considered in ref.[5].) 

In fact, the expression (2.1) is generated as the effective action in QED in curved 
spacetime including vacuum polarisation effect s[l]. It is valid in the approximation of 
weak curvature and low frequency photons, so that the neglect of higher powers of the cur- 
vature tensor and extra covariant derivatives in the interaction terms is justified. Roughly, 
(see ref.[13] for a careful discussion), this effective action is a good approximation in the 
parameter range A c < A < L, where A is the photon wavelength, L is a typical curvature 
scale, and A c = 1/m is the electron Compton wavelength, m being its mass. The electron 
provides the appropriate quantum scale as the effect is generated by vacuum polarisa- 
tion diagrams involving an internal electron loop. In this case, the coupling constants are 
given to one loop in terms of the fine structure constant by a = — j^i^, b = 3^— and 
c = — 3gQ^- Notice, however, that if this is the origin of the action (2.1), then the va- 
lidity of the results we deduce here concerning photon propagation is only established for 
relatively low frequency photons, whereas to consider the speed of real signal propagation 
and the implications for causality we need to consider high frequencies. The question of 
dispersion, the generalisation of the effective action to include higher derivative terms and 
the physical interpretation of the results will be discussed in ref.[13]. 

The Bianchi identity and equation of motion following from the action (2.1) are* 

D^F U \ + D U F X ^ + DxFp,, = (2.2) 

and 

D„F^ - -^(^bR^xD^ + AcR^ Xp D^F xp ^j = (2.3) 

where we have assumed that both the photon wavelength and m _1 are small compared to 
the curvature scale. 

The simplest way to determine the characteristics of photon propagation from these 
equations is to use geometric optics. In the leading geometric optics approximation, we 
write the electromagnetic field strength as the product of a slowly varying amplitude and 
a rapidly varying phase, i.e. F pv = f^ u exp iO, where the wave vector is k^ = 3^9. In the 
quantum particle interpretation, we identify as the photon momentum. The amplitude 
is constrained by the Bianchi identity to be of the form f pu = k^a^ — k v a p , where the 
direction of a M specifies the polarisation. Light rays (photon trajectories) are defined as 
the integral curves of the wave vector (photon momentum). Without the additional term 
in eq.(2.3), these may easily be shown to be null geodesies. This is no longer true when 
the equivalence principle violating interactions are included, nor is the light cone condition 
k 2 = necessarily satisfied. 

Our main concern is with the local modifications to the light cone induced by these 
direct curvature interactions. Introducing an orthonormal frame using vierbeins defined 
by e%e „, where r\ ab is the Minkowski metric, the photon equation of motion 

(2.3) gives 

k 2 a b - ^(^R ac {k a k c a b - k a k b a c ^ - ^ ^R abcd k a k c a d ^ = (2.4) 

* See ref. [12] for an extended version of this section. 
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The polarisation vector is spacelike normalised, a b cib = — 1, and can be taken orthogonal 
to the momentum, k b ab = 0. Of the three remaining degrees of freedom, in the quantum 
theory only the two transverse polarisation states are physical. Given a polarisation vector 
a b satisfying this equation, the corresponding light cone condition is simply 

k 2 - ^-R ac k a k c + ^-R ahcd k a k c a b a d = (2.5) 
We shall now explore the consequences of these equations for the Kerr geometry. 



3. The Kerr Spacetime 

The Kerr spacetime (see, e.g., refs. [14-16]) is described by the metric 

rfs 2 = _ p 2 + p2 + p2dQ2 + J_ E 2 s[n 2 e S d(j) _ udt y (31) 

£ 2 A p 2 

where 

2aMr 9 9 9 9 „ 

<jj = ^ 2 p z = r 2 + a 2 cos 2 9 

A = r 2 -2Mr + a 2 E 2 = (r 2 + a 2 ) 2 - a 2 sin 2 9 A (3.2) 

This metric represents the exterior spacetime of a rotating black hole. It is axially sym- 
metric about the rotation axis 9 = 0. It is specified by two parameters, M and a, where 
M is the mass and Ma the angular momentum as measured from infinity. For a = 0, it 
reduces to the Schwarzschild spacetime. In the description below, we assume a < M. 

The condition A(r) = 0, for which there is a coordinate singularity similar to that 
in the Schwarzschild metric, has two solutions, r = r± = M ± \JM 2 — a 2 . The larger, 
r = r + , is in fact the event horizon, from within which no particle on a timelike or null 
trajectory can escape to infinity. The region r < r_ contains a ring singularity. Just as in 
the Schwarzschild geometry, there is a Killing vector which is timelike in the asymptotic 
region (large r) but which is spacelike within the event horizon. However, for the Kerr 
spacetime, there is a further region outside the horizon, r+ < r < te(9), for which the 
Killing vector remains spacelike. This is the ergosphere. te{9) = M + \/M 2 — a 2 cos 2 9 

is the value of r for which the metric component g tt = — ^1 — ^f 2 ^ vanishes and the 

Killing vector is null. The outer limit of the ergosphere, the stationary limit surface, is 
the inner boundary of the region where particles travelling on timelike curves can remain 
at rest relative to infinity. Within the ergosphere, even null curves are pulled round in 
the direction of the rotation. Particles may, however, escape to infinity from this region. 
The situation is summarised pictorially in Fig.l, taken from ref.[15]. The stationary limit 
surface and event horizon coincide at the poles, where the stationary limit surface becomes 
null. At the equator, is equal to 2M, the Schwarzschild radius. 

We now introduce a local orthonormal frame. The appropriate basis 1-forms are e a 
(a = 0, 1,2,3) with 



e 



= e° t dt e 1 = e 1 r dr e 2 = e 2 9 d9 e 3 = e% (dcf) - udt) (3.3) 
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where the vierbeins are 



e° t = -p^ e\ = -^= e 2 e = p e 3 t = -we 3 ^ = — ^Esin# (3.4) 



The Kerr metric is Ricci flat, so R a b = 0. There are six independent non-vanishing 
components of the Riemann curvature which we can choose in this frame to be[16] 

-Roioi = A -R0202 = B -R0123 = C 

-R0231 = D -R0102 = E -R0113 = F (3-5) 

The complete set of non- vanishing components is 

-R2323 = — -R0101 -R1313 = — -R0202 -R1212 = — -R0303 = -R0101 + -R0202 

-R0312 = — -R0123 — -R0231 -R3132 = -R0102 -R0223 = — -R0113 (3.6) 

together with those related by the usual symmetries of R a bcd- The quantities A, . . . , F are 
given by 



A = ^(r 2 -3a 2 cos 2 e)^(2(r 2 +a 2 ) 2 + a 2 Asm 2 e) 
(r 2 - 3a 2 cos 2 9) ± ((r 2 + a 2 ) 2 + 2a 2 A sin 2 e) 



Mr 



A + B = ^L(r 2 -3a 2 cos 2 9) 



aM cos 9 
C = - 

aM cos 9 



(3r 2 - a 2 cos 2 9) ± (2(r 2 + a 2 ) 2 + a 2 A sin 2 9 
(3r 2 - a 2 cos 2 9) ± ((r 2 + a 2 ) 2 + 2a 2 A sin 2 9) 



=> C + D = - aM ™ Se (3r 2 -a 2 cos 2 9) 

E = - aM J S 9 (3r 2 - a 2 cos 2 9) ^ (r 2 + a 2 ) sin 9 

TP Mr (n 2 2 2 M 3fl V^/ 2 , 2\ ■ /> (1 T\ 

F = — g- (3r — a cos 0) — — j — (r + a )sm0 (3.7) 

Introducing the notation = 5®Sl — S^S^ etc., we can rewrite the complete Riemann 
tensor compactly in the following form: 

R abcd = 2A (y%u£ - u 2 Ju 2 J - u$u£ + u%u%) 

+ 2b(u 02 U 02 - U 13 U 13 - U 03 U 03 + U l2 U 12 ^) 

+ c(u 01 u 23 + u 23 u 01 - u 03 u 12 - U l2 U 03 ) 
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+ d(-u 02 u 13 - U 13 U 02 - U 03 U 12 - U 12 U 03 ) 

+ e(u 01 u 02 + u 02 u 01 + u 13 u 23 + u 23 u 13 ) 

+ f(u 01 U 13 + U 13 U 01 - U 02 U 23 - U 23 U 02 ) (3.8) 

where we have suppressed the a, 6, c, d indices after the first line for clarity. To obtain 
the curvature components Rp V \ p in the coordinate frame (£, r, 6, 4>) , simply replace by 
U®1 = e° At e 1 !y — e 1 I/ e° At etc. in this expression. 

Notice the following simplifications for special cases. On the event horizon, A(r+) = 0, 
so E = F = 0, A = —2B and C = —2D. In the equatorial plane 6 = ir/2, cos6> vanishes 
and so C = D = E = 0. 



4. Photon Propagation in Kerr Spacetime 

We now return to eq.(2.4) describing photon propagation in curved spacetime. For 
Ricci flat spacetimes such as Kerr, this reduces to 

k 2 a b + eR abcd k a k c a d = (4.1) 

where we have written e = — 8c/m 2 . This is a set of three simultaneous linear equations 
for the independent components of the polarisation a&. To solve these [1], it is convenient 
to introduce the following linear combinations of momentum components: 

l b = k a U% m b = k a U° a 2 n b = k a U°J (4.2) 

together with the dependent combinations 



Pb = 


k a Ul 2 = 


-^{k l m b - 


- k 2 l b ) 




Qb = 


k a Ul 3 = 


■j-o( kl n b - 


k 3 i b ) 




r b = 


k a u 2 J = 


^o(k 2 n b - 


k 3 m b ) 


(4.3) 



The vectors £, m, n are independent and orthogonal to k a . 

We can therefore rewrite (4.1) as a set of equations for the independent polarisation 
components a.£, a.m and a.n by contracting appropriately. Substituting eq.(3.8) for the 
Riemann tensor, we therefore arrive at the following set of equations: 

= k 2 a.£ + 2Ae[f a.£ - £.r a.r - £.n a.n + £.p a.p^j 

+ 2Be^£.m a.m — £.q a.q — £.n a.n + £.p a.p^j 
+ Ce(^£ 2 a.r + £.r a.£ — £.n a.p — £.p a.nj 
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+ De(^—i.m a.q — i.q a.m — i.n a.p — i.p a.n^j 
+ Ee(i 2 a.m + i.m a.i + i.q a.r + i.r a.q^j 
+ Fe[i 2 a.q + i.q a.i - i.m a.r - i.r a.rnj 



2 

= k a.m + 

= k 2 a.n + ... (4.4) 

To save space, we have only written the first in full. The second two have a similar form 
and are easily reconstructed. 

In principle, this set of equations could now be solved in general. However, it is much 
more illuminating to look at a selection of special cases which illustrate the most important 
features. First, consider photon propagation in the equatorial plane. 



(i) Equatorial plane, radial motion 

To illustrate the method of solution, consider first radial motion confined to the 
equatorial plane, where C = D = E = 0. The photon momentum components satisfy 
k 2 = k 3 = 0. The various momentum-dependent terms appearing in eqs.(4.4) therefore 
simplify considerably. In particular, we have i 2 = k°k° — fe 1 fc 1 , m 2 = n 2 = k°k° and 
i.m = i.n = m.n = 0, while for the others we find p 2 = q 2 = k*k l , m.p = n.q = k°k x , with 
all other contractions vanishing. For the polarisation projections, a.p = k 1 /k°a.m, a.q = 
k 1 /k°a.n and a.r = 0. Substituting these special results into eqs.(4.4) and rewriting the 
system in matrix form, we find 

'k 2 +2Ati 2 F ^ 2 w 

k 2 + 2Aek 1 k 1 + 2Be(k°k° + k 1 k 1 ) 

Fe^k 1 k 2 - 2Aek°k° - 2Be(k°k° + k 1 ^) 




x a.m = (4.5) 



In general we would have to diagonalise to find the polarisation eigenvectors, the 
corresponding values of k 2 being given as solutions of the vanishing of the determinant. 
In this case, however, there is a further simplification. We should regard e as a small 
parameter (in fact, it is the dimensionless combinations Ae, Be etc. which are small) since 
the original action involving single powers of the curvature will be valid only for A c < L, 
and work to consistent order in small e. Since with e = the light cone condition is just 
k 2 = 0, in general we have k 2 = 0(e). Now, for radial motion i 2 = —k 2 so is of O(e), 
and therefore the off-diagonal entry proportional to F is actually of 0(e 2 ) and should be 
neglected at lowest order. 

We therefore find the solutions: 

k 2 = 0, corresponding to the polarisation a a proportional to i a , 
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k 2 + 2Aek°k° + 2Be(k°k° + k l k l ), corresponding to a a proportional to m a , 

k 2 — 2Aek°k° — 2Be(k°k° + k l k l ), corresponding to a a proportional to n a . 

The solution with k 2 = corresponds to an unphysical polarisation, l a = — /c 1 ^, while 
the two physical transverse polarisations are proportional to m a = k°5 2 , i.e. polarisation 
in the 9 direction, and n a = k°5 3 , i.e. polarisation in the <fi direction. At this level the 
effect is non-dispersive and the (phase or group) velocity shift is simply 



Sv 



k ( 



-1 = ±(A + 2B)e = T e(^^A) (4.6) 



k 

for 9 ((p) polarisation respectively. 

Two features of this result are immediately apparent. First, the shift in k 2 , or equiva- 
lently the velocity shifts away from the conventional speed of light, are equal and opposite 
for the two transverse polarisations. Second, the shift vanishes on the event horizon, since 
A(r+) = 0. Both these observations turn out to be examples of general theorems and we 
discuss them further in the next section. 



(ii) Equatorial plane, orbital motion 

Now consider photons travelling in the orbital ((f)) direction in the equatorial plane. 
In this case, k 1 = k 2 = 0. The analysis goes through in the same way as described above 
and gives the solutions: 

k 2 + 2Aek°k° - 2Bek 3 k 3 - 2Fek°k 3 = 0, corresponding to a a proportional to £ a , 

k 2 - 2Aek 3 k 3 + 2Bek°k° + 2Fek°k 3 = 0, corresponding to a a proportional to m a , 

k 2 — (2 A + 2B)ek°k° = 0, corresponding to a a proportional to n a . 

The two physical transverse polarisations, in the radial and 9 directions, are given 
by the first two of these solutions respectively. Notice, however, that now the light cone 
condition is not quadratic in each of the momentum components separately due to the 
presence of the F term. This introduces a splitting in the velocities for propagation with 
k 3 > and k 3 < 0. For the radial polarisation we find 

Sv = (A-B T F)e = (^(r 2 + a 2 ) 2 + a 2 A T Mr 2 + a 2 )VA) (4.7) 

depending on whether the motion is with or against the direction of spin. For the 9 
polarisation, Sv has the opposite sign. This expression simplifies on the horizon, where we 
have 

X I 3M I A 

Hhorizon = —3- e ( 4 - 8 ) 

r + 

and on the stationary limit surface = 2M, where (for k 3 < 0) 

Hstat. lim. = ^4(^2^) ( r i + 5a2 ( r l + a2 )) 6 ( 49 ) 
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(iii) Radial motion, arbitrary direction 



In the Schwarzschild and Reissner-Nordstrom spacetimes, the light cone for radial 
photons remains k 2 = for any direction (the solutions are of course spherically symmetric) 
and for any value of the radial coordinate r. In the Kerr spacetime, however, we have just 
seen that k 2 7^ for radial photons in the equatorial plane except at the horizon r = r + . 
We now want to check whether this result remains true independent of the angle 9 to the 
polar axis at which the photons are directed. 

The calculation follows the lines of (i) except that now the curvature components C, 
D and E are non- vanishing. Eq.(4.5) generalises to: 

k 2 + 2Ae£ 2 Ee£ 2 Fe£ 2 ^ \ 

Eek°k° k 2 + 2Aek 1 k 1 + 2Be(k°k° + k 1 k 1 ) -(C + 2D)ek°k 1 

Fek^k 1 -{C + 2D)ek°k 1 k 2 - 2Aek°k° - 2Be(k°k° + k x k x ) J 

(al\ 

x a.m =0 (4.10) 
y a.n J 

Again there is a solution k 2 = corresponding to the unphysical polarisation a a propor- 
tional to l a . The transverse polarisations diagonalising (4.10) are complicated r and 9 
dependent linear combinations of the unit vectors in the 9 and <p directions, the corre- 
sponding velocity shifts being 

Sv = ±( K (A + 2B) 2 + l -(C + 2D) 2 Y 

= ±eA 3MaWfl ^ 2(r2 _ 3fl2 cQg2 Q)2 + 1 fl2 cQg2 ^ (3r2 _ fl2 cQg2 Q)2 ^ l (4n) 

Notice that 5v = along the polar axis, while of course the previous result is recovered in 
the equatorial plane. Most important, however, we again see that the velocity shifts are 
equal and opposite for the transverse polarisations and that, independently of 9, Sv = 
on the event horizon where A(r_|_) vanishes. 



5. The Polarisation Sum Rule, the Horizon and the Stationary Limit Surface 

These results on modifications of the light cone in special cases, together with those 
previously obtained in refs[l,2], motivated the formulation of two general theorems which 
were stated precisely and proved in ref.[12] (see also ref.[17]). These are: 

Polarisation Sum Rule 

This relates the sum over the transverse polarisation states to an appropriate projection 
of the Ricci tensor, viz. 

Vfc 2 = — \r(Ab + Sc)R ac k a k c (5.1) 
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Horizon Theorem 



This states that at the event horizon, the light cone for photons with momentum directed 
normal to the horizon remains k 2 = 0, independent of the polarisation. 



Clearly, the specialisation of the polarisation sum rule to Ricci flat spacetimes, i.e. 
^ pol k 2 = 0, or equivalently X) po i ^ v = 0> * s satisfied by all the examples in Kerr spacetime 
discussed in section 4. 

Similarly, for radial photons, eq.(4.12) ensures that k 2 = 0, or equivalently v = 1, 
at the event horizon even though the velocity differs from 1 for all other values of the 
radial coordinate. The horizon theorem therefore ensures that the geometric event horizon 
remains a true horizon for photons propagating according to the action (2.1), even in a 
spacetime with as rich a structure as Kerr. 

In the light of this, it is interesting to ask whether the geometric stationary limit sur- 
face r = te{0) specified by gtt = retains its defining property for real photon propagation. 
Recall that this is the surface within which even light signals emitted against the direction 
of rotation are pulled round, through the phenomenon of dragging of inertial frames[14,15], 
so that as measured by an asymptotic observer they propagate in the direction of rotation 
of the black hole. 

The results of section 4(ii) show that this is not true. The fact that the light cone is 
modified for photons emitted in the negative direction (see eqs.(4.7)(4.9)) even on the 
stationary limit surface shows that the effective stationary limit surface is shifted from 
t\e(6>) to a larger or smaller value of r depending on the photon polarisation. 

The shift is most readily found using a trick used in ref. [1] to calculate the modification 
to the bending of light in a Schwarzschild spacetime. As we have seen, the light cone 
condition is modified to fj a i,k a k b = where fj a b differs from the orthonormal (Minkowski) 
metric by the terms of 0(e) calculated in section 4. Propagation with this modified light 
cone in a spacetime with metric g^ v is therefore equivalent to conventional propagation 
with light cone rj a bk a k b = in a modified geometry with metric g^ v = fj a be a ^e b u . 

In this case, we have determined fj a fj for propagation in the counter-orbital direction 
in the equatorial plane (section 4(ii)). Since g u = —e t e t + e 3 t e 3 t, we see that evaluated 
on the (unperturbed) stationary limit surface e° t = e s t . We therefore have immediately 
that 



9tt 



g tt + 2e(A -B + F)e° t e° t 



(5.2) 



fl=f,r=2M 



Setting g tt = gives the effective stationary limit surface. In the equatorial plane, we 
therefore find 



te 



eff 



2M±e(A-B + F) 



a 2 + 2M 2 



(5.3) 



depending on the polarisation, the magnitude of the velocity shift term being given in 
eq.(4.9). 

Finally, we should emphasise again that these results follow from taking the action 
(2.1) literally. If instead we regard it as the lowest-order effective action for QED, we should 
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ask to what extent the inherent approximations can be relaxed. The weak curvature 
(A c < L) approximation is difficult to improve on, so the magnitude of the results will 
necessarily be extremely small in the domain where the derivation is reliable. This limits 
the interest of these results for astrophysical black holes. This is, however, just the usual 
situation for quantum field phenomena in curved spacetime such as, e.g., the Hawking 
effect [18]. The predicted effects are of the order of a quantum scale divided by a curvature 
scale, so are expected to become large only for microscopic black holes or in the very early 
universe. A more serious limitation of the action (2.1) is the implicit restriction to relatively 
low-frequency photons (A c < A), which poses serious questions as to the observability and 
relevance for signal propagation and causality of the shifts in the light cone considered here. 
These issues will be addressed in a forthcoming paper [13]. Despite the need for caution, 
however, the general theorems inferred from our study of photon propagation in black 
hole spacetimes may well prove to be valid outside these limitations. In particular, the 
horizon theorem, whose proof relied on the physical result that classically no gravitational 
radiation (or matter) crosses the horizon[19], looks sufficiently robust to conjecture that it 
may be true in general. 
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Fig. 1 A picture of the equatorial plane of the Kerr geometry with a 2 < M 2 . The circles 
represent the position after a short time interval of photons emitted in the r, <p plane from 
the points represented by the heavy dots. The distinctive 'frame-dragging' effect giving 
rise to the ergosphere is readily seen, as is the property of the horizon that it cannot be 
crossed even by null trajectories. 
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